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We studied the nonequilibrium magnetisation reversal in kinetic Ising ferromagnets driven 
by periodic impulsive magnetic field in meanfield approximation. The meanfield differ- 
ential equation was solved by sixth order Runge-Kutta-Felberg method. The periodicity 
and strength of the applied impulsive magnetic field play the key role in magnetisation 
reversal. We studied the minimum strength of impulsive field required for magnetisation 
reversal at any temperature as a function of the periodicity of the impulsive field. In the 
'^ ■ high temperature and small period this is observed to be linear. The results are compared 

with that obtained from Monte Carlo simulation of a three dimensional Ising ferromagnet. 
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I. Introduction: 

The Ising model is a prototype to study some equilibrium as well as nonequilibrium 
phenomena [1] . Recently, the nonequilibrium aspects of Ising ferromagnets in presence of 
time varying fields become an interesting field of modern research [1]. The time dependent 
meanfield equation [2], of kinetic Ising ferromagnet in presence of applied magnetic field, 
is a simple representation of nonequilibrium phenomena. The nonequilibrium dynamic 
phase transition [1] is studied extensively in this system. It exhibits tricritical behaviour 
and a method of finding the tricritical points is proposed recently [3]. Like so called 
equilibrium transition, the nonequilibrium dynamic transition shows the divergence of 
time scale and power law divergence was found recently [4]. 

The magnetisation reversal, by nucleation, in the Ising ferromagnet was studied by 
Monte Carlo simulation (for large systems) in the presence of constant (in time) magnetic 
field [5]. This magnetisation reversal is an equilibrium phenomenon. Due to the presence 
of thermal fluctuation, any small amount of field (constant in time) may lead to the mag- 
netisation reversal. In that sense no such minimum amount of field is required. However, 
in the meanfield approximation a minimum amount of field is required for magnetisation 
reversal. The externally applied periodic (in time) impulsive fields keeps the system al- 
ways away from equilibrium. So, it would be interesting to know how the nonequilibrium 
magnetisation reversal occurs if one applies a periodic impulsive fields. In this paper, we 
addressed this question and studied the nonequilibrium magnetisation reversal (in kinetic 
Ising ferromagnets in meanfield approximation) in presence of periodic (in time) impul- 
sive magnetic fields. The nonequilibrium magnetisation reversal was studied [6] recently 
in Ising model both by Monte Carlo simulation and by solving the dynamical meanfield 
equation. The dynamics of magnetisation reversal was also studied [7] experimentally in 
polycrystalline Co film. 

We have arranged the paper as follows: the meanfield differential equation of kinetic 
Ising ferromagnet and the method of numerical solution are described in the next section. 
The numerical results are given in section-Ill. The paper end with a summary in section- 
IV. 



II. Model and numerical solution: 

The differential equation of instantaneous average magnetisation m{t) of kinetic Ising 
ferromagnet, driven by a time varying magnetic field, in meanfield approximation, is given 

as [2] 

dm , ,m + h(t). ,_,, 

T^ = -m + tanh( ^), (1) 

where, h{t) is the externally applied time varying magnetic field and T is the temperature 



measured in units of the Boltzmann constant (Kb)- This equation describes the nonequi- 
hbrium behaviour of instantaneous value of magnetisation m{t) of Ising ferromagnet in 
meanfield approximation. Here, r stands for the microscopic relaxation time for the spin 

flip [2]. 

The periodic impulsive magnetic field is given as 



h{t) = — /io,for t = nAt,where n = 1, 2, 3. 

= otherwise. (2) 

Where At and Hq denote the periodicity and strength of the impulsive field respectively. 
We have solved this equation by sixth order Runge-Kutta-Felberg (RKF) [8] method 
to get the instantaneous value of magnetisation m{t) at any finite temperature T, Hq and 
At. The method of solving the ordinary differential equation ^ = F(t,'m(t)), by sixth 
order RKF method, is described briefly as: 

m(f J- riA — m(i\ -^ fl^fci , 6656fc3 , 28561^4 _ 9^5 i 2^6^ 
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ki = dt- F{t,m{t)) 



k2 = dt 
ks = dt 
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The time interval dt was measured in units of r (the time taken to flip a single spin). 
Actually, we have used dt = 0.01 (setting r=1.0). The local error involved [8] in the sixth 
order RKF method is of the order of {dt)^{= 10^^^). We started with initial condition 
m(t = 0) = 1.0. 



III. Results: 

Starting from the initial magnetisation mlt = 0) = 1.0, the instantaneous magnetisa- 
tion (m(t)) decreases as time goes on. At any instant the magnetisation m{t) becomes 
negative and this is known as magnetisation reversal. We studied this magnetisation 
reversal by applying the periodic impulsive fleld (Eqn.-2). The time required for mag- 
netisation reversal depend on the value of At and ho and the temperature T and is called 
reversal time T^. All times are measured in the unit of O.Olr throughout this study. A 
typical results of magnetisation reversal is shown in Fig.-l, for At = 10 and Hq = 1.5 
at temperature T = 0.90. Here, the magnetisation reversal occurs at time t = 760. It 



was observed that, at any temperature T, the reversal occurs early as we increase the 
strength ho keeping At fixed. In this case, a minimum value of strength, h^^^, is required 
for reversal. We checked few cases, taking the value of ho less than h'^^^ and observed no 
reversal even for t = 5 x 10^. 

At any fixed temperature T, the minimum field strength (/ij„j„) required for mag- 
netisation reversal, decreases as the periodicity (At) decreases. Similarly, for fixed At, 
/i^j„ decreases as the temperature increases. In the present study, we observed that, at 
temperature T = 0.80, for At = 6 and At = 8 the values of /i^j„ becomes 0.33 and 0.43 
respectively. Similarly, at temperature T = 0.90, for At = 6 and At = 8 the values of 
^min become 0.13 and 0.17 respectively. It may be noted here that for static (in time) 
field the minimum values calculated for equilibrium magnetisation reversal are 0.07 and 
0.03 for T = 0.80 and T = 0.90 respectively. At any particular temperature, the value 
of minimum field required (for periodic impulsive fields in nonequilibrium case) is always 
higher than that for the static field. These results are shown in Fig.-2. In this figure, 
the reversal time Tj. is plotted against the field strength ho for different values of At and 
temperature T. The reversal time Tr is observed to increase here with the decrease of 
field strength ho- 

At any fixed temperature T, the minimum required field strength /ij„j„, increases with 
the periodicity At of the impulsive fields. In the high temperature (T^ = 1 for equilibrium 
ferro-para transition) limit, it is observed that the relationship is linear (for small value 
of At). We have shown this in Fig. -3. 

These results are compared with that obtained from Monte Carlo simulation. We 
have considered a three dimensional Ising model with nearest neighbour ferromagnetic 
interaction under periodic boundary condition in all directions. The Hamiltonian of this 
system is represented as 

H = —JT,^ijySiSj — h{t)lliSi (3) 

where Sj = ±1 is Ising spins, J is the nearest neighbour ferromagnetic interaction strength 
and h{t) is periodic impulsive field (Eqn-2). We have considered a cubic lattice of size 
L = 20 (here). The initial condition was taken as all spins are up (i.e.. Si = +1 for all 
i). The updating scheme is taken as random. We select a spin randomly and calculated 
the energy in fiipping, i.e., ^E. We fiipped the selected spin with probability P/ = 
Min[l,exp(— Aii^/i^sT")], following the Metropolis algorithm[9]. Here Kb is Boltzmann 
constant. The temperature T is measured in the unit of J / Kb- L^ number of such 
random updating of Ising spins constitutes a single Monte Carlo step. This defines the 
unit of time in this simulation. In three dimensions the ferro-para transition temperature 
is Tc = 4.511 [9]. Now we allowed the system to follow the Metropolis dynamics in the 
presence of periodic impulsive field h(t)- At any fixed temperature T, we have calculated 
(averaged over 10 different random updating sequences) the minimum field required h"^^^ 



for the reversal of magnetisation and studied it as a function of the interval or periodicity 
At of periodic impulsive field. Here also, we have observed the linear variation of /ij„m 
with At. We have studied this for two different temperatures, here, T = 4.0 and T = 4.2. 
We observed that as the temperature increases the slope of the straight line (/i^j„ versus 
At) decreases. This is shown in Fig-4. These Monte Carlo study supports the meanfield 
results. 

IV. Summary: 

We have studied the nonequilihrium magnetisation reversal of kinetic Ising ferromag- 
nets in presence of periodic impulsive magnetic fields in meanfield approximation. At any 
finite temperature and fixed value of the periodicity of the impulsive field a minimum 
value of field strength is required for magnetisation reversal (sign change of instantaneous 
magnetisation). This value of minimum required field strength depends on the periodicity 
of the impulsive fields. As the periodicity increases the minimum field required increases. 
In the high temperature limit, it is observed to be linear in the periodicity (for small 
values) of the applied impulsive fields. This study differs from the case where the mag- 
netic field is constant in time and consequently this leads to equilibrium magnetisation 
reversal. However, in this case since the field is time dependent, the characteristics of 
the magnetisation reversal is nonequilibrium type. These are studied in Ising meanfield 
dynamics. As a comparison the the nonequilibrium magnetisation reversal is also studied 
in three dimensional Ising ferromagnet driven by periodic impulsive fields by Monte Carlo 
simulation. The Monte Carlo results and the meanfield results are in good agreement. 
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Fig-1. The time variation of magnetisation m(t) and the apphed periodic impulsive field 
h{t) of periodicity At = 10 and strength ho = 1.5. The magnetisation reversal happens 
at t = 760. Here, the temperature T = 0.90. 
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Fig-2. The variation of magnetisation reversal time (T^.) with respect to the strength of 
impulsive field (/iq) for different temperatures (T) and periodicity (At) of the impulsive 
fields. 
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Fig-3. The variation of /i]^j„ with At for different temperatures (T) represented by 
different symbols. T = 0.75(0), T = 0.80(n), T = 0.85(A) and T = 0.90(0). Continuous 
straight hues represent the hnear best fit. 
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Fig. -4. The plots of /ij„j„ versus At for two different temperatures, obtained from Monte 
Carlo simulation. The errorbars in each data are shown by small vertical lines. The 
continuous straight line represents the linear best fit. 
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